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SUMMARY 

Trajectories of water droplets moving in the ideal two-dimensional 
flow field ahead of a body of rectangular cross section and infinite ex- 
tent in the downstream direction have been calculated by means of a dif- 
ferential analyzer. Data on collection efficiency and distribution of 
water impingement are presented. 


INTRODUCTION 

As part of a comprehensive research program dealing with the prob- 
lem of icing protection for aircraft, an investigation of the impingement 
of cloud droplets on airfoils and other aircraft components has been un- 
dertaken by the NACA Lewis laboratory. In addition to specific aircraft 
components (refs. 1 to 3), certain general aerodynamic bodies such as 
cylinders (ref. 4) and ellipsoids of revolution (ref. 5) have been studied 
because various aircraft components may be approximately represented by 
•these shapes. Among the general body shapes of interest is that of a 
rectangular body having a flat surface facing the air stream. Such a body 
may be used as an approximation to the shape of the sensing elements of 
certain instruments used to measure icing-cloud properties. 

Cloud-droplet size and size distribution are sometimes measured by 
collecting samples of cloud droplets on small transparent slides coated 
with oil or some water-sensitive substance (refs. 6 to 8) . The slides 
are then photographed to obtain a record of the droplet size distribution. 
Various types of sampling devices have been used in which the collecting 
element is usually a flat surface perpendicular to the airstream. The 
evaluation of the cloud-droplet size distribution from samples obtained 
in this way requires a knowledge of the relation between droplet diameter 
and collection efficiency. DrOplet-implngement data applicable to a body 
of rectangular cross section would also be useful in estimating the col- 
lection efficiency of the central portion of the rotating-disk icing-rate 
meter (refs. 7 to 9) and the rectangular bars used to measure the spray 
distribution in an icing wind tunnel (ref . lO) . 
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Trajectory calculations based on ideal two-dimensional flow about a 
ribbon (ref. ll), which is essentially a rectangle of zero chord, have 
been used as an approximation to obtain the impingement on rectangular 
bodies of finite chord. However, the ideal flow field ahead of a ribbon 
is not a very satisfactory representation of the actual flow field in 
front of a rectangular body because of the dependence of the flow ahead 
of the body on conditions existing behind the body. 

The ideal flow field about a ribbon is symmetrical fore and aft, as 
shown in figure l(a), with the streamlines closing in directly behind the 
body. In the case of real flow about a rectangular body, on the other 
hand, the streamlines cannot be calculated exactly; but general knowledge 
of the flow about solid boundaries suggests that the flow field is some- 
tMng like that shown schematically in figure l(b) . Flow separation oc- 
curs at the edge, and a wake is formed extending downstream. The effect 
of variations in the ratio of chord to thickness cannot be estimated ac- 
curately, but is probably small, since the wake is roughly equivalent to 
an indefinite downstream extension of the body. 

Potential flow exists outside the limits of the boundary layer and 
wake (shown approximately by the dashed line). Figure l(c) shows the 
ideal-flow streamlines around a rectangular half body (a body of infinite 
extent downstream). This flow field resembles the estimated real flow in 
figure 1(b) mudi more closely than does the ideal flow about a ribbon 
(fig. 1(a)). The effect of the downstream extent of the body on the flow 
field ahead is shown in figure l(d), which is a superposition of portions 
of the fields shown in figures l(a) and (c) . The comparison shows that 
the presence of the afterbody exerts a considerable influence on the flow 
field ahead of the body. 

The ideal flow field of the rectangular half body represents a closer 
approximation to real flow than does the ideal flow about a ribbon; how- 
ever, it too is subject to the disadvantage that the velocity becomes in- 
finite at the corner. In spite of this disadvantage, the ideal flow field 
in front of the rectangular half body was used, in the calculation of tra- 
jectories, because it was considered to be sufficiently realistic to yield 
approximate results of practical value and because the velocity components 
could be calculated without too much difficulty. 


SYMBOLS 

The following symbols are used in this report: 
A area per unit span, sq ft/ft 

a droplet radius, ft 

B,C empirical constants, dimensionless 
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Cjj coefficient of drag, dimensionless 

E collection efficiency, dimensionless 

f functional relation between Yq and yg 

G flux of liquid water per unit span, slugs/ (sec) (ft) 

K inertia parameter, dimensionless 

L half width of rectangle, ft 

n impingement distribution index, dimensionless 

p,q arbitreiry numbers 

Re local Reynolds number with respect to droplet, dimensionless 

Re^ free-stream Reynolds number with respect to droplet, dimensionless 

t time , sec 

U free-stream velocity, ft/sec 

Ujj.,Uy local air velocity components, ratio to free-stream velocity, 
dimensionless 

Vjj,Vy local droplet velocity components, ratio to free-stream velocity, 
dimensionless 

W * total rate of water Impingement per unit span, slugs/ (sec) (ft) 

local rate of water impingement per unit area, slugs/ ( sec ) (sq ft) 

U 

w free-stream liquid-water content, slugs/cu ft 

x,y rectangular coordinates, ratio to half width L, dimensionless 

Y distance from centerline (Y = Ly), ft 

y^ trajectory starting ordinate, dimensionless 

y trajectory ordinate at point of impingement, dimensionless 

s 

3 local impingement efficiency, dimensionless 
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3o local impingement efficiency at centerline 

3-]^ local impingement efficiency at edge 

ji viscosity of air^ slugs/(ft) (sec) 

density of air, slugs/cu ft 

density of water, slugs/cu ft 

T dimensionless time function, t = ^ 

XJ 

2 

® Re^/K, dimensionless 


ANALYSIS 

In order to find the rate and distribution of droplet impingement 
on the surface of a body, it is necessary to determine the cloud-droplet 
trajectories with respect to the body. The method used in this case to 
calculate the trajectories of cloud droplets is described fully in ref- 
erence 4. Assumptions that are necessary to the solution of the problem 
are: (l) At a large distance ahead of the body the droplets are at rest 

with respect to the air; (s) the only external force acting on the drop- 
lets is the drag force due to the relative velocity of the air with re- 
spect to the droplets; and (3) the droplets are always spherical and do 
not change in size. 


Differential Equations of Droplet Motion 

The differential equations that describe the motion of the droplets 
were obtained by equating the drag force with the rate of change of mo- 
mentum of the droplet (ref. 4). The equations are expressed in dimension- 
less form in order to maintain the number of calculations at a minimum 
and to simplify the presentation of the results. They apply to the mo- 
tion of droplets in a plane perpendicular to the edges of the rectangular 
half body in the system of rectangular coordinates shown in figure 2. 

The origin is located at the center of the front face of the rectangle, 
and the unit of distance is the half width L of the rectangle. The 
dimensionless coordinates x and y are ratios to L, The dimension- 
less air velocity components u and u and droplet velocity components 

X y 

v^ and Vy are ratios to the free-stream velocity U, which is the unit 

of velocity. Time is expressed in terms of the dimensionless quantity 
T = tU/L. The unit of time is the time required to go a distance L at 



NACA TN 3658 


5 


a speed U. At a large distance ahead of the body, a uniform air flow 
bearing cloud droplets is assumed to be approaching the rectangle from 
the negative x direction moving parallel to the x axis at a speed U 
(fig. 2) . 


The equations of motion of the droplets are 


dT 



— K K - ’'x) 




V 


^ i (u - V ) 

24 K '■ y y'' y 


(1) 


where Re is the local Reynolds number with respect to the droplet diame 
ter 2a and the local relative velocity between the air and the droplet, 
thus 


Re . - ''x'^ 

The coefficient in equation (2) is called the free-stream Reynolds number 
ReQ) hence , 


Reo = 


2apg^U 


( 3 ) 


The dimensionless number, K, the inertia parameter, is defined as follows 

(4) 


9 pL 


The coefficient of drag Cj^ may be obtained from experimental data as a 
function of the local droplet Reynolds number Re. 


An examination of equations (l) and (2) shows that the characteris- 
tics of the trajectories depend only on the values of the two dimension- 
less parameters K and Rsq. Thus, a unique set of droplet trajectories 

exists for every combination of K and Rcq. 


Method of Solution 

The differential equations of motion (l) are difficult to solve be- 
cause values of the velocity components and the factor containing the 
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coefficient of drag depend on the position and velocity of the droplet at 
each instant and, therefore, are not known until the trajectory is traced. 
The values of these quantities must be fed. into the equations as a tra- 
jectory is developed. This was accomplished by using a mechanical dif- 
ferential analyzer constructed at the MCA Lewis laboratory for this pur- 
pose (ref . 4) . The results were obtained in the form of plots of droplet 
trajectories in the coordinates of figure 2. 


The following information was required for use with the differential 
analyzer in the trajectory calculation: 

(1) Coefficient of drag data: Values of the factor Cj^ Re/24 as a 

function of Re were obtained from tables in reference 11. 

(2) Air. velocity components: The air velocity components u^ and 

Uy were determined as functions of x and y by means of an analytical 

solution of the two-dimensional potential flow field ahead of a rectangu- 
lar half body. Because of symmetry with respect to the centerline (see 
fig. 1), one half of this field is identical to the flow field around a 
step discontinuity in the floor of a channel of Infinite depth, a configu- 
ration discussed in reference 12. equations used to calculate the 

velocity components were obtained by means of the Schwarz-Christoffel 
transformation as described in reference 12. These equations, expressed 
in parametric form, are as follows: 




- P cosh P - sinh p cos q 

sin^q + sinh^p 

cosh p sin q - sin q cos q 

" 2 2 
sin q + sinh p 


> 


X = zi (p + sinh p cos q) 


( 5 ) 


n 


y = i (q + cosh p sin q) 

rt 




Any pair of values of the arbitrary parameters p and q deter- 
mines the coordinates x and y of a point in the flow field and the 

components u and n of the velocity at that point. Values of u . 

X y ^ x’ 

Uy, X, and y were calculated from equation (5) for a number of points 

in the flow field ahead of the rectangle, and curves giving u^ and Uy 

as functions of x and y were obtained by cross-plotting. The calcu- 
lated values of coordinates and air velocity components are presented in 
table I. 
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(3) Starting conditions: The integrations vere started at x = -60 

where the air velocity con^jonents are = 0.9945 and Uy = O.OOOly 

for all values of y between 0 and 1. The initial droplet velocity 
components were assumed to be = 0.9945 and Vy = 0. The calcula- 
tions were performed in three steps: x = -60 to -lOj x = -10 to -Ij 

and X = -1 to 0. Different scale factors were used on the machine for 
each step, and the results of each of the first two steps were used to 
establish starting conditions for the succeeding step. 


Determination of Droplet Impingement 

The differential analyzer traces the trajectories of the droplets 
from a starting point defined by the ordinate at a large distance 

ahead of the body to the final point of impingement having the ordinate 

y at the surface of the rectangle. The trajectories provide pairs of 
s 

values of Yq and y^ that establish the relation between yQ and yg 

for various combinations of K and RCq. The relation is expressed sym- 
bolically as follows: 


Vq = f(ys) 


( 6 ) 


and the symbols f (O) and f (l) are used to denote the values assumed by 

y^ when y^ is 0 and 1, respectively, 
u s 


If a rectangular body of half width L moves through a cloud of 
uniform liquid-water content and uniform droplet size, the flux of liquid 
water per unit span between any two trajectories with starting ordinates 
Ly^ ^ and lyQ 2 (^yQ dimensions of length) is 


' 0,1 


G = UwL(yo^2 ' yo,l) 


Since the trajectories define the path of the flow of liquid water, the 
same flux arrives at the rectangle surface between the trajectory terminal 
ordinates Ly -1 and Ly „ corresponding to the starting ordinates 

X 2 

Ly^ T and Ly« The surface area per unit span of the strip of the 
rectangle between the two terminal ordinates Ly_ and Ly _ is 

S ^ X S ^ o 

A = i'(ys,2 - 

The liquid-water flux G is distributed over the area A} hence, the 
average rate of water impingement per unit area on this strip is 



8 


NACA TN 3658 


G 

A 


UwL(y, 

L(y 


0,2 


s,2 


' ^ 0,1 
^ 3 , 1 ^ 


) 


If is allowed to approach yg ^2 as a limit, 

"by G/a is the local rate of water impingement W 

y 

Y = y_L from the centerline. 

o 


the limit approached 
at the distance 




= lira ^ = \Jw 
A-^0^ 


% 

<^3 


The local impingement efficiency 3 is defined as follows: 


(V) 


Hence , 



3 




( 8 ) 


( 9 ) 


The total rate of water interception per unit span over the entire 
face of the rectangle is twice the integral of the local impingement rate 
from the stagnation point (Y = O) to the outside edge (Y = L) ; thus. 



W = 2LUw[f (1) - f (0)] 

where the symbol f refers to the function of equation (6) . 


In a symmetrical flow field, a droplet approaching the stagnation 
point along the central streamline does not change direction; thus, the 
value of yQ corresponding to y^ = 0 is 0; therefore. 


W = 2LUwf (l) 


( 10 ) 


The collection efficiency E 


is defined as follows: 


W 

2LUw 


( 11 ) 
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Hence , 

E = f(l) (12) 

The collection efficiency is therefore equal to the value of Yq when 
yg = 1, which is the nondimensional starting ordinate of the trajectory 

that strikes the outside edge of the rectangle in the coordinate system 
of figure 2. 


RESULTS 

A total of 74 pairi of values of yQ and yg were calculated with 
the differential analyzer using 23 combinations of Roq and K. It was 
found that the relation of equation (6) between yQ and yg for any 
particular combination of K and RCq could be represented approximate- 
ly by an empirical equation of the fonn 

yo = ®ys + cys ' 

where B and C are positive and less than 1 and n is positive and ' 
greater than 1. The collection efficiency E is the value of Yq when 

yg = 1; hence. 


E = B + C (14) 

The local impingement efficiency ^ may be obtained by differenti- 
ating equation (13) as follows: 




(15) 


The local in 5 >ingement efficiency at the stagnation point Pq is obtained 
by setting yg = 0 in equation (15); thus,' 


and from equation (14), 


C = E - 

Substituting in equation (13) gives 


yQ = + (E - 3o)ys 


(16) 
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and from equation (15) there is obtained 

P = pQ + n(E - Po)y|'^'^^ (17) 

Equation (17) shows that the local impingement efficiency has a minimum 
value p^ at the stagnation point and increases to a maximum given by 

Pi = Po + n(E - Po) (18) 


at the outside edges of the rectangle. 

The Impingement pattern corresponding to any given pair of values 
of K and RSq is determined by three peireimeters : (l) the collection 

efficiency E, (2) the local impingement efficiency at the stagnation 
point Pq, and (3) the distribution index n, which determines the form 

of the distribution of impingement. 


An analysis of the trajectory data has shown that these three param- 
eters, which are functions of K and Re^, may be represented approxi- 
mately by the following empirical equations: 


-logePo = (j 


JL 

n = 1 + e K + [l.l + (loggO.72 K)^] 

0.538\A 


(19) 


K + 10 


+ 0.2 + 0.0238 Re> 


87 


/ \ 115 + ROq + 0.004 Reg 


) 


(K - 0.28)"°‘® 


( 20 ) 


-log E = 


151 


K'" + 150 


+ 0.267 + 0.225 Re 




180 

250 + R6 




(K - 0.15)"°-’^^ (21) 

Equations (17), (19), (20), and (2l) provide a means of determining 
approximately the collection efficiency and the distribution of impinge- 
ment for any combination of Rcq and K. The trajectory calculations 

used in deriving the equations adequately covered a range of K from 
1/2 to 10 and of RSq from 16 to 256. Outside this range the equa- 
tions represent extrapolations. 
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Values of n, and E calculated from the empirical equations 

are presented in figures 3, 4, and 5, respectively. The 74 values of 
yQ obtedned with the differential analyzer and the corresponding values 

calculated from the equations are presented for comparison in table II. 

The standard deviation of the differences between corresponding values 
is 0.0051; therefore, the probable error involved in using the equations 
to determine yQ is ±0.0034 based on the assumption that the results from 

the trajectories are absolutely correct. The dimensionless parameters 
RSq and K are defined in equations (3) and (4).. Equations and graphs 

for use in determining values of K and Roq in terms of air speed, al- 
titude, body size, and droplet diameter in practical units are presented 
in reference 13. 


DISCUSSION 

The use of empirical equations to represent the impingement patterns 
is essentially a means of intei^iolating and extrapolating the trajectory 
data. The procedure is analogous to that of drawing a family of curves 
and interpolating between them. Because of random errors in the trajec- 
tory data, a certain amount of smoothing is necessary in using either 
procedure. The use of empirical equations has the suivantage of accomplish- 
ing the smoothing operation in a consistent manner over the entire range 
of K and ROq, thus assuring that the results are as accurate for inter- 
mediate values of these parameters as for the values for which trajecto- 
ries were calculated. The equations also provide a more reliable basis 
for extrapolation to values of K and Rsq beyond the range of the 
trajectory calculations . 

The average accuracy of the equations is probably a little better 
than is indicated by the probable error of ±0.0034 in yQ based on the 

departures in table II, because these deviations are due partly to errors 
in the trajectory calculations. The same degree of accuracy (probable 
error about 0.003) also applies to values of E and from equations 

(20) and (2l). Values of p at intermediate points and values of p-^ 

at the corner of the rectangle, calculated from equations (17) and (18), 
may be subject to leirger errors because of the effect of errors in the 
distribution index n. 

The validity of the representation of the impingement distribution 
patterns given by equation (17), using values of n, P^, and E from 

equations (19), (20), and (21), respectively, can be checked by a com- 
parison with average values of p for the intervals of yg between 

adjacent trajectories. Such values are obtained directly from the tra- 
jectory data by means of the relation 
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Ay, 

Ay 


( 22 ) 


T1i6 comparison is shown in figure 6 in which the smooth curves were 
calculated from equation ( 17 ), and the block diagrams represent averages 
from equation (22)' based on individual intervals between adjacent cal- 
culated trajectories. An examination of figure 6 reveals the general 
pattern of the impingement distribution as a function of K and Rcq, 

At high values of K, the impingement efficiency approaches 1 aoross the 
entire face of the rectangle, but more slowly for high values of RCq. 

At low values of K, the impingement efficiency approaches 0, but the 
number of trajectories calculated (3 at K = l/3) is too small to deter- 
mine accurately the shape of the impingement profile. For all values of 
K, the impingement efficiency decreases, and the distribution becomes 
more uniform with increasing ReQ. The greatest curvature of the impinge- 
ment profile, indicating a sharp maximum impingement at the edge, occurs 
at low values of Re^ and values of K between 1 and 2. The curves cal- 
culated from equation (17) agree satisfactorily with the values of 3 ob- 
tained from adjacent trajectories over the entire range of K and Roq 

for which trajectories were determined. Equation (17), therefore, pro- 
vides a consistent representation of tte variations of the impingement 
profile with K and 

It is of interest to compare the droplet-impingement data for a 
rectangle with similar data for a ribbon in reference 11. The results 
are presented in reference 11 in the form of curves showing E and 3 q 
as functions of K and $ ($ = Re^/K) . Values of E and 3^ read 
from these curves are compared in table III with values calculated from 
equations (16) and (l7) for corresponding values of K and Roq. For 

most values of K and $ the rectangle has a lower collection effi- 
ciency than the ribbon^ this effect is due to the influence of the after- 
body in reducing the rate of deceleration of the air in the region just 
ahead of the body. 

It was pointed out in the discussion of flow fields in the INTRODUC- 
TION that the flow field of the rectangular half body differs from that 
of the ribbon by the addition of an afterbody, but that both differ from 
the real flow about a rectangular body because of the absence of flow 
separation and because of the requirement of infinite velocity at the 
edges. An examination of the comparative streamlines in figure 1(d) shows 
a great difference in velocity in the immediate vicinity of the edge. 
(Velocity is inversely proportional to streamline spacing.) In spite of 
the difference in the flow fields, however, the differences in collection 
efficiency of the ribbon and rectangle are only moderate. The maximum 
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difference in table IH is 0.099 at K = 2 and 0 = 10,000, and the 
average difference for all values in table III is 0,026. These facts 
suggest that,the results obtained using the ideal flow ahead of the rec- 
tangular half body are likely to be a fairly good approximation to the 
actual impingement on a rectangular body in real flow. The largest dif- 
ference between the ideal flow and the real flow is that the ideal flow 
permits a high local air velocity in the neighborhood of the edge; how- 
ever, these differences in local air velocity in the neighborhood of the 
edge have a relatively small effect on the trajectories. 


CONCLUDING REMAEKS 

The rate and distribution of water impingement on a semi-infinite 
rectangular body moving in a cloud of uniform liquid-water content and 
droplet size have been determined by means of droplet trajectory calcula- 
tions. Empirical equations have been found providing an approximate rep- 
resentation of the impingement rate and distribution in terms of the 
inertia parameter K and the free-stream droplet Reynolds number ReQ. 

Equation (l7) gives the local impingement efficiency p as a function of 
distance from the centerline of the rectangle in terms of the following 
three parameters: (l) the collection efficiency E, (2) the impingement 

efficiency at the centerline p^, and (3) the impingement distribution 

index n. 


These parameters may be determined from K and Rcq by means of 
empirical equations (19), (20), and (2l) . 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 14, 1955 
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TABLE I. - VELOCITY FIELD AHEAD OF RECTANGULAR HALF BODY 

fThe table gives values of the dimensionless space 
coordinates x and y and corresponding values 
of dimensionless velocity components u^ and 

u^ (see fig. 2) calculated from equation '(5) ^ 


X 

y 


u 

y 

X 

y 



-0.09567 

0 


0.07485 

0 


0 

0.12551 

0 

0.09918 

-.19242 



. .14889 



-.09473 

.12622 

.07571 

.09855 

-.29136 



. 22128 



-.19053 

.12835 

.15036 

.09699 

-.39364 



.29132 



- . 28847 

. 13195 

. 22335 

.09429 

-.50049 



.35837 



-.38969 

.13712 

.29394 

.09073 

-.61323 



.42191 



-.49538 

. 14395 

.3615 

.08637 

-.73335 



.48154 



-.60686 

.15260 

.4253 

.08142 

-.86246 



.53707 



-.72556 

. 16328 

.4852 

.07603 

-1.0518 



.60437 



-.85310 

.17622 

.5408 

.07040 

-1.2655 



.66404 



-1.0400 

.19750 

.6081 

.06275 

-1.6451 



.73977 



-1.2507 

. 22419 

.6677 

.05523 

■ -2.1190 



.80050 



-1.6248 

.27674 

.7432 

.04469 

-2.9586 



.86173 



-2.0914 

. 34869 

.8033 

.03542 

-4.1437 



.90515 



-2.9171 

.48638 

.8639 

.02536 

-6.3798 



.94138 



-4.0816 

.69271 

.9067 

.01779 

-9.9599 



.96403 



-6.2771 

1.0996 

.9424. 

.01119 

-11.96 



.97045 



-9.7905 

1.7711 

.9647 

.00695 

-14.40 



.97580 



0 

. 24681 

0 

.19892 

-17.30 



.98010 



-.09202 

. 24818 

.07782 

.19774 

-20.9 



.98368 



- . 18505 

.25233 

.15466 

.19434 

-25.2 



.98660 



- . 28009 

.25935 

. 22960 

.18880 

-30.5 



.98903 



-.37821 

. 26940 

.30184 

.18144 

-37.0 



.99100 



-.48056 

.28271 

.3707 

.17251 

-44.8 



.99264 



-.58836 

. 29957 

.4356 

.16237 

-54.5 



..99397 



-.70299 

.32036 

.4960 

.15138 

-66.0 



.99505 



-.82589 

.34556 

.55204 

.13995 

-98 



.99668 



-1.00567 

.38700 

.61937 

.12446 

-145 



.99777 



-1.2079 

.43897 

.67846 

.10929 

-216 

% 

r 

.99851 



-1.5659 

.54132 

.75277 

.08814 

-320 

1 


.99900 



-2.0110 

. 68141 

.81161 

.06968 
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TABLE I . - Continued. VELOCITY FIELD AHEAD OF RECTANGULAR HALF BODY 

[The table gives values of the dimensionless sp^e 
coordinates x and y and corresponding values 
of dimensionless velocity components and 

Uy (see fig. 2) calculated from equation (5) .] 


X 

y 


u 

y 

X 

y 


. ^ 

-2.7964 

0.94955 

0.87027 

0.04975 

-0.14923 

0.58990 

0.19037 

0.52027 

-3.9010 

1.3514 

.91143 

.03482 

-.22535 

.60516 

.25295 


0 

.36350 

0 

.30257 

-.30333 

.62702 

.36594 

.47837 

-.08763 

.36549 

.08318 

.30063 

-.38384 

.65597 

.44476 

.45015 

-.17617 

.37150 

.16220 

. 29528 

-.46762 

769261 

^.51652 

.41876 

-.26652 

.38167^ 

.24046 

. 28647 . 

_-. 55549 

.73784 

.58097 

.38561 

-.35965 

.39623 

.31554 

. 27476 

-.64834 

-.79265 

..-.63827 

.35190 

-.45658 

.41550 

.38666 

. 26064 

-.78167 

.88280 

.7'6392 

.30766 

-.55842 

.43993 

.45316 

. 24473 

-.92848 

.99580 

■ .75851 

. 26574 

-.66641 

.47005 

.51469 

.22756 

0 

;-7 66908 

0 

.66820 

-.78186 

.50656 

.57037 

. 20947 

-.06609 

.67239^ 

. ' 

^'0801 

- .66394 

-.95012 

.56659 

.63835 

.18583 

-.13259 

.68241 

.21326 

.64699 

-1.1386 

.64187 

.69669 

.16258 

-.19992 

.69936 

31323 

.62186 

-1.4706 

.79014 

.76902 

.13044 

-.26854 

.7236L 

'.40600 

.58914 

-1.8804 

.99310 

.82538 

.10266 

-.33890 

.75573 

,.49028 

.55077 

-2.6012 

1.3815 

.88073 

.07293 

-.41152 

.79644 

.56537 

--^50883 

0 

.47508 

0 

.41422 

-.48696 

.84664 

. 63121 

.46506 

-.08163 

.47762 

.08763 

.41151 

-.56585 

.90750 

.68822 

.42123 

- . 16403 

.48529 

.17377 

.40350 

-.67760 

1.0075 

.75165 

.36466 

- . 24798 

.-49826 

.25708 

.39066 

0 

.74914 

0 

.81949 

-.33428 

.51683 

.33641 

.37364 

-.05716 

.75266 

.12466 

.81181 

-.42382 

.54141 

.41084 

.35325 

-.11454 

.76331 

.24522 

.78960 

-.51753 

.57258 

.47965 

.33043 

-.17235 

. 78129 

.35810 

.75456 

-.61646 

.61100 

.54257 

.30599 

-.23081 

.80704 

.46069 

.70953 

-.70243 

.65756 

.59954 

.28086 

- . 29017 

.84113 

.55148 

.65756 

-.87424 

.73412 

.66631 

. 24740 

-.35069 

.88433 

.62994 

.60170 

-1.0440 

.83015 

.72327 

: 21528 

-.41266 

.93758 

.69645 

.54458 

-1.3404 

1.0193 

.79232 

.17142 

-.47639 

1.0022 

.75206 

.48848 

0 

.57788 

0 

.53542 

0 

.81831 

0 

1.0000 

-.07431 

.58087 

.09618 

.53158 

-.04775 

.82190 

.14888 

.98885 
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lEABLE I. - Concluded. VELOCITY FIELD AHEAD OF RECTANGULAR HALF BODY 

[The table gives values of the dimensionless space 
coordinates x and y and corresponding values 
of dimensionless velocity components u^^ and 

liy (see fig. 2) calculated from equation (5) .] 


X 

y 



X 

y 



-0.09549 

0.83274 

0.29131 

0.95663 

-0.02118 

0.95495 

0.32951 

1.8215 

-.14324 

.85109 

.42190 

.90665 

-.04176 

.96398 

.62023 

1.6952 

-.19099 

.87735 

.53705 

.84356 

-.06113 

.97924 

.76417 

1.3668 

-.23873 

.91212 

.63516 

.77240 

-.07864 

1.0011 

1.0088 

1.3188 

-.28648 

.95617 

.71632 

.69781 

0 

.97508 

0 

2.4142 

- .33423 

1.0105 

.78179 

.62349 

-.01386 

.97762 

.49501 

2.3248 

0 

.87506 

0 

1.2203 

-.02696 

.98529 

.90034 

2.0906 

-.03833 

.87858 

.18480 

1.2034 

- .03850 

.99826 

1.1755 

1.7863 

-.07645 

.88923 

.35876 

1.1552 

0 

.98942 

0 

3.3050 

-.11413 

.90721 

.51336 

1.0817 

-.00786 

.99141 

.85669 

3.1543 

-.15116 

.93296 

.64377 

.99151 

-.01482 

.99742 

1.4292 

2.6018 

-.18729 

.96705 

.74884 

.89289 

0 

.99681 

0 

5.0275 

- . 22227 

1.0103 

.83013 

.79296 

-.00347 

.99818 

1.7229 

4.3791 

0 

.91908 

0 

1.4966 

0 

.99959 

0 

10.078 

-.02940 

.92239 

.23954 

1.4687 

-.00076 

1.0003 

4.8931 

6.3865 

-.05840 

.93241 

.45956 

1.3943 





-.08656 

.94936 

.64606 

1^2826 





-.11343 

.97361 

.79306 

1.1508 





-.13856 

1.0057 

.90171 

1.0130 





0 

.95196 

0 

1.8675 
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TABLE II. - COMPAEISON OF CALCULATED STARTENG ORDINATES 

[starting ordinates of trajectories calculated with differential 

analyzer and corresponding values of Yq calculated from empiri- 
cal equations for same values of yg are compared.] 


Reo K 



yo 

(Trajectory) 


^ yo 

Difference 

(Equations) 


0.227 

-0.001 

.461 

0 

0.037 

+0.004 

.070 

+ .008 

0.074 

-0.001 

.149 

0 

.225 

+ .001 

.257 

-.008 

0.082 

0 

.166 

+ .003 

.250 

+ .005 

.327 

0 

.405 

-.003 

.495 

+ .003 

.503 

-.002 

0.266 

+0.007 

.355 

+ .009 

.528 

+ .005 

.697 

-.005 

0.186 

0 

.371 

+ .002 

.557 

0 

.745 

+ .005 

0.384 

0 

.916 

0 

0.012 

-0.001 

0.037 

0 

.062 

+ .008 

.095 

+ .007 

.124 

+ .015 


0.676 

1.028 


0.228 

.461 



1 

0-182 

0.082 


.367 

.163 


.552 

.245 


.715 

.327 


.856 

.408 


.986 

.492 


.996 

.505 

2‘ 

0.414 

0.259 


.551 

.346 


.807 

.523 


1.012 

.702 

4 

0.238 

0.186 


.472 

.369 


.703 

.557 


.926 

.740 
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TABLE II. - Continued. COMPAEISON OF CALCUIATED STARTING ORDIMTES 

[Starting ordinates Yq of trajectories calculated with differential 
analyzer and corresponding values of Yq calculated from empiri- 
cal equations for seime values of Yg are compared.] 


Reo 

K 


yo 

(Trajectory) 

yo 

(Equations) 

Difference 

64 

1 

0.127 

0.040 

0.037 

-0.003 



.490 

.158 

.146 

-.012 



.722 

.237 

.229 

-.008 



.974 

.359 

.352 

-.007 



1.035 

.400 

.390 

-.010 


2 

0.500 

0.251 

0.254 

40.003 



.644 

.338 

.333 

-.005 



.920 

.511 

.511 

0 


4 

0.261 

0.180 

0.178 

-0.002 



.519 

.361 

.361 

0 



.765 

.542 

.545 

+.003 



.985 

.727 

.720 

-.007 


10 

0.443 

0.378 

0.378 

0 



.977 

.857 

.856 

-.001 

128 

1 

0.5 

0.545 

0.035 

0.034 

-0.001 



.820 

.053 

.060 

+.007 


1 

0.184 

0.038 

0.037 

-0.001 



.351 

.077 

.072 

-.005 



.515 

.114 

.110 

-.004 



.638 

.153 

.143 

-.010 



.910 

.232 

.233 

+ .001 


2 

0.399 

0.164 

0.164 

0 



.576 

.248 

.244 

-.004 



.748 

.331 

.331 

0 



.931 

.440 

.436 

-.004 


4 

0.287 

0.176 

0.174 

-0.002 



.565 

.352 

.354 

+ .002 



.700 

.440 

.446 

+ .006 



.944 

.621 

.621 

0 
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TABLE II. - Concluded. CCMPARISON OF CALCULATED STARTING ORDINATES 

[starting ordinates of trajectories calculated vith differential 

analyzer and corresponding values of Yq calculated from empiri- 
cal equations for same values of yg are compared.] 


Reo 

K 

^s 

yo 

(Trajectory) 

yo 

(Equations) 

Difference 

128 

10 

0.461 

0.372 

0.373 

+0.001 



.903 

.751 

.753 

+.002 

256 

0.5 

0.832 

0.033 

0.027 

-0.006 


1 

0.266 

0.035 

0.036 

40.001 



.517 

.072 

.075 

+ .003 



.730 

.110 

.115 

+.005 



.943 

.148 

.161 

+.013 


2 

0.483 

0.159 

0.160 

+0.001 



.697 

.243 

.243 

0 



.975 

.367 

.366 

-.001 


4 

0.466 

0.259 

0.251 

-0.008 



.634 

.343 

.349 

+ .006 



.905 

.520 

.516 

-.004 


10 

0.488 

0.369 

0.370 

40.001 



.947 

.746 

.745 

-.001 




MCA TN 3658 


21 


TABLE III . - COMPARISON OF COLLECTION EFFICIENCY E AND STAGMTION- 
LINE IMPINGEMENT EFFICIENCY ’ FOR RECTANGLE AND CORRESPONDING 


DATA FOR RIBBON FROM REFERENCE 11 (LANGMUIR) 


K 

6 ) 

E 

(Rectangle) 

E 

(RLbbon) 

(Rectangle) 

(Ribbon) 

0.25 

0 

0.123 

0.122 

0 

0 


100 

.060 

.084 


0 


i;ooo 

.036 

.062 

0 

0 


10,000 

.011 

.035 

0 

0 

0.333 

0 

. 0.262 

0.244 

0.052 

0.098 


100 

.161 

.168 

.034 

.076 


1,000 

.113 

.125 

.017 

.063 


10,000 

.049 

.075 

.003 

.046 

0.50 

0 

• 0.436 

0.450 

0.287 

0.235 


100 

.316 

.350 

.231 

.186 


1,000 

.244 

.260 

.161 

,154 


10,000 

.134 

.157 

.066 

.114 

1 

0 

0.651 

0.666 

0.554 

0.470 


100 

.535 

.595 

.484 

.380 


1,000 

.452 

.503 

.384 

.320 


10,000 

.305 

.382 

.228 

.244 

2 

0 

0.789 

0.808 

0.722 

0.667 


100 

.691 

.734 , 

. .650 

.574 


1,000 

.610 

.666 

.548 

.492 


10,000 

.466 

.565 

.388 

.389 

4 

0 

0.878 

0.894 

0.835 

0.800 


100 

.802 

.836 

.769 

.731 


1,000 

.733 

.783 

.676 

.660 


10,000 

.610 

.692 

.539 

.544 

10 

0 

0.953 

0.954 

0.928 

0.909 


100 

.903 

.918 

.874 

.864 


1,000 

.855 

.883 

.805 

.823 


10,000 

.773 

.819 

.709 

.737 

20 

0 ' 

0.982 

0.976 

0.962 

0.951 


100 

.948 

.954 

.923 

.921 


1,000 

.915, 

.929 

.872 

.890 


10,000 

.862 

.882 

.803 

.834 

100 

0 

0.997 

0.995 

0.992 

0.989 


100 

.982 

.988 

.970 

.978 


1,000 

.968 

.980 

.946 

.965 


10,000 

.948 

.960 

.912 

.942 





Figure 1. - Comparison of two-dimensional flow fields. 





Figure 2. - Coordinate system and notation for droplet trajectory calculations. 
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c 



.2 .4 .6 1 2 4 6 10 20 40 60 100 

Inertia parameter, K 


Figure 3. - Impingement distribution index calculated from equation (19) ♦ 




Collection efficiency 
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Inertia parameter, K 

Figure 5. - Collection efficiency calculated from equation (21). 
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